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1. Introduction 

The Bertrand-Darboux theorem is well known to provide a necessary and sufficient 
condition for two-degree-of-freedom natural Hamiltonian systems associated with the 
Euclidean metric to admit an integral of motion quadratic in momenta (Darboux 
1901). Moreover, that condition is necessary and sufficient for the natural Hamiltonian 
systems to be separable in either Cartesian, polar, parabolic or elliptic coordinates 
(Marshall and Wojciechowski 1988). This theorem thereby provide a sufficient 
condition for the complete integrability of the natural Hamiltonian systems, which 
will be referred to as the Bertrand-Darboux integrability condition (BDIC) in the 
present paper. The BDIC has been studied repeatedly from various viewpoints; the 
separation of variables (Marshall and Wojciechowski 1988, Grosche et al 1995 in the 
path-integral formulation), the complete integrability (Perelomov 1990), and so-called 
the direct method (Hietarinta 1987), for example (see also Whittaker 1944 as an older 
reference and the references in the above-cited literatures). 

The aim of the present paper is to show that the BDIC is obtained as an outcome 
of a study on the Birkhoff-Gustavson(BG)-normalization (Gustavson 1966, Moser 
1968) for the Hamiltonians of the perturbed harmonic oscillators with homogeneous 
CM&zc-polynomial potentials (PHOCP's). In solving an inverse problem of the BG- 
normalization for a given PHOCP-Hamiltonian with a help of computer algebra, the 
family of the perturbed harmonic oscillators with homogeneous gttarfic-polynomial 
potentials (PHOQP's) is identified which share the same BG- normal form up to degree- 
4 with the given PHOCP. Consequently, a new deep relation is found between the 
BDIC for the PHOCP's and that for the PHOQP's: It is shown that if a PHOCP 
and a PHOQP share the same BG-normal form up to degree-4 then both oscillators 
are integrable in the sense that they satisfy the BDIC. It is worth noting that the 
present work was inspired in the debugging process of the computer-programme named 
'ANFER' for the BG- normalization (Uwano et al 1999, Uwano 2000), where the one- 
parameter Hcnon-Heiles system was taken as an example. 

Before the outline of the present paper, the inverse problem of the BG- 
normalization is explained very briefly. The BG-normalization has been a powerful 
method for nonlinear Hamiltonian systems. For example, when a two-degree-of- 
freedom Hamiltonian system with a l:l-resonant equilibrium point is given, the 
BG-normalization of its Hamiltonian around the equilibrium point provides an 
'approximate' Hamiltonian system: The truncation of the normalized Hamiltonian 
up to a finite degree is associated with the approximate Hamiltonian system, which 
provide a good account of the surface of section with sufficiently small energies 
(Kummer 1976, Cushman 1982). Such a good approximation implies that finding 
the class of Hamiltonian systems admitting the same BG-normalization up to a finite- 
degree amounts to finding a class of Hamiltonian systems which admit the surface of 
section similar to each other. The following question has been hence posed by the 
author as an inverse problem of the BG-normalization (Chekanov et al 1998, 2000, 
Uwano et al 1999): What kind of polynomial Hamiltonians can be brought into a given 
polynomial Hamiltonian in BG-normal form ? Since elementary algebraic operations, 
differentiation, and integration of polynomials have to be repeated many times to solve 
the inverse problem, computer algebra is worth applying to solve the inverse problem; 
see Uwano et al (1999) and Uwano (2000) for the programme named 'ANFER' and 
Chekanov et al (1998, 2000) for 'GITA -1 '. 

The organization of this paper is outlined as follows. Section 2 sets up the 
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ordinary and the inverse problems of the BG-normalization for Hamiltonians, which 
will be often referred to as 'the ordinary problem for Hamiltonians' and 'the inverse 
problem for BG-normal form Hamiltonians', respectively, henceforth. Although the 
setting-up of the ordinary problem seems to be merely a review of Moser (1968), it 
is of great use to define the ordinary problem in an mathematically-sound form : 
Through the review, the class of canonical transformation to be utilized in the BG- 
normalization can be specified explicitly. In section 3, the one-parameter Henon-Heiles 
Hamiltonian is taken as an example to illustrate how the ordinary and the inverse 
problem are proceeded. From the solution of the inverse problem, it follows that if 
the one-parameter Henon-Heiles system and the perturbed harmonic oscillator with 
a homogeneous-quartic polynomial potential (PHOQP) share the same BG-normal 
form up to degree-4 then both dynamical systems satisfy the BDIC, so that they 
are integrable. In section 4, the discussion in section 3 made for the one-parameter 
Henon-Heiles Hamiltonian is extended to the Hamiltonians of the perturbed harmonic 
oscillators with homogeneous-cubic polynomial potentials (PHOCP's): The ordinary 
and the inverse problems of the BG-normalization for the PHOCP Hamiltonians 
are dealt with there. In the ordinary problem, the Hamiltonians for PHOCP's are 
normalized up to degree-4. For the BG-normal form thus obtained, the inverse problem 
is solved up to degree-4. Both of the problems are solved with (a prototype of) the 
symbolic-computing program ANFER working on Reduce 3.6 (Uwano 2000). It is 
shown that if a PHOCP and a PHOQP share the same BG-normal form up to degree- 
4 then both oscillators are integrable in the sense that they satisfy the BDIC. Section 5 
is for concluding remarks including a conjecture on a further extension. 

2. Setting-up the ordinary and the inverse problems of 
BG-normalization 

In this section, the ordinary and the inverse problems of the BG-normalization are 
reviewed for the two-degree-of-freedom Hamiltonians (Uwano et al 1999). 

2.1. The ordinary problem 

We start with defining the ordinary problem along with a review of Moser (1968). Let 
R 2 x R 2 be the phase space with the canonical coordinates (q,p) (q,p 6 R 2 ), and 
K(q,p) be the Hamiltonian expressed in power-series form, 



where each Kk(q,p) (k = 3, 4, • • ■) is a homogeneous polynomial of degree-fc in (q,p). 

Remark 1. The convergent radius of the power series (|l]) may vanish. This happens to 
any K that is not analytic but differentiable around the origin, for example. In such 
a case, the power series (Q) is considered only in a formal sense. We will, however, 
often eliminate the word 'formal' from such formal power series henceforth. 

We normalize the Hamiltonian K through the local canonical transformation of the 
following form. Let (£, 77) be another canonical coordinates working around the origin 
of R 2 x R 2 . We consider the canonical transformation of (q,p) to (£,ri) associated 




(1) 



with a generating function (Goldstein 1950) in power-series form, 



W{q,r,)=Y,<lim + ^W k {q,ri), (2) 

j=l fc=3 

where each W k is a homogeneous polynomial of degree-fc in (q, rf) (fc = 3, 4, ■ ■ •). The 
W(q, 17) is said to be of the second-type since W(q, n) is a function of the 'old' position 
variables q and the 'new' momentum ones r\ (see Goldstein 1950). The canonical 
transformation associated with W(q, rf) is given by the relation 

(<7,PJ -► ISjf?) Wltn P= ~dq ^~drf' ^' 

which leaves the origin of R 2 x R 2 invariant on account of (|2|). By G(£,rf), we denote 
the Hamiltonian brought from K{q 1 p) by 



G Wv =K v^) (4) 

through the transformation (^|). The BG-normalization for K is accomplished by 
choosing the generating function W(q,rf) in (Q) to put G(£, 77) in BG-normal form: 

Definition 2.1. Let G{fc,rf) be written in the power-series form, 

2 00 

G(t, r,) = - + & + £ G k & V ), (5) 

i'=l k=3 

where each G k (t;,r]) is a homogeneous polynomial of degree-fc (k = 3, 4, • • •) in (£, 77). 
The power series G(£, 77) is said to be in BG-normal form up to degree-r if and only if 



(6) 



holds true, where {•, -}^. n is the canonical Poisson bracket (Arnold 1980) in (£, 77). 

Let us equate the homogeneous-polynomial part of degree-fc (fc = 3, 4, • • •) in (jij) . Then 
equation (Q) is put into the series of equations, 

G k (q,r 1 ) + (D q , n W k ) = K k (q,r 1 )+3> k (q,r ] ) (fc = 3,4,---), (7) 

where D q ^ is the differential operator, 



The $ k (q,rf) ln (|3) is the homogeneous polynomial of degree-fc in [q, rf) which 
is uniquely determined by W3, ■ ■ ■ , Wk-i, Ks,---,K k ^i, G 3 , ■ ■ ■ , Gfe_i given: In 
particular, we have $3(5,77) = and 



^ 1 (dW 3 y dK 3 



7 = 1 



dW- 
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Since we will deal with only the BG-normalization up to degree-4 in the present paper, 
we will not get the expression of <J>fc for k > 4 into more detail (see Uwano et al 1999, 
if necessary). 

To solve equation (Q), the direct-sum decomposition induced by D qv of the 
spaces of homogeneous polynomials is of great use. Let Vk(q,r)) denote the vector 
space of homogeneous polynomials of degree-fc in (q, rj) with real-valued coefficients 
(fc = 0, 1, • • •). Since the differential operator D qJ1 acts linearly on each Vk(q,r]), the 
action of D q ^ naturally induces the direct-sum decomposition, 

V k (q, rj)= imaged W©ker£>« (k = 0, 1, ■ • •), (10) 
of T4 (q, 77) , where Dq^ denotes the restriction, 

D ^ = D i,v\v klq ,n) ( fc = 3 > 4 '"-)- (U) 

(k) (k) 

Remark 2. We have kerD^ = {0} and imageDg^ = Vk(q, 77) if k is odd. 

For D qv and (k = 3, 4, • ■ •), we have the following easy to prove. 
Lemma 2.2. Equation (JsJ) is equivalent to 

(D q . v (G k I (,,,)))(«, r,) = (D (*> (G fc I fej7) ))( 9 , 77) - (fc = 3, • • • , r). (12) 

Namely, Gfc|( 9l ,) G ker-D^ (fc = 3, • • • , r). 

According to (|l(]), let us decompose K)~{q, 77) and "3?/c(g, 77) (fc = 3, 4, • • •) to be 
K k (q, r?) = »?) + V), 

Kr*\q, 77), S^fo V) G imaged, 
^ er (g,7 ? ),$r(?,'7)GkerDW,. 



where 



Since Gfe G kerDg^ by Lemma 2.2 and since D q ^Wk G imagel?^ , we obtain 

G fe ( g ,7 ? ) = ^ er ( g ,r7) + ^ er (g,r7) (fc = 3,4,---), (15) 
as a solution of (|t]), where is chosen to be 

W k {q,rj) = (d^ ^ (Kr s %,r,)+^r%,r,))) (16) 

What is crucial of ( |l6| ) is that Wfc G LmagaD^ (fc = 3,4, • • •): For a certain integer 
K > 3, let us consider the sum, = W K + (any polynomial in kerD^i), where G k 
and Wk with k < n are given by (|l5|,[l6|). Even after such a modification, 1F K satisfy (^) 
with k = k still, which leads another series of solutions of (^) with k > k. Therefore, 
under the restriction, Wk G imageZ?^ (fc = 3, 4, • • ■), we can say that ( [l5| ) with ( |l6| ) is 
the unique solution of (0). To summarize, the ordinary problem is defined as follows: 
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Definition 2.3 (The ordinary problem). For a given Hamiltonian K(q,p) in power 
series (H), bring K{q,p) into the BG-normal form G(£,77) in power series (^J) which 
satisfy (Q) and (^) with r = oo, where the second- type generating function W of the 
form (||) is chosen to satisfy ([|) and 

W k (q, r?) e imaged $ (fc = 3,4, • • •). (17) 

Theorem 2.4- The BG-normal form G(£, 77) for the Hamiltonian K(q,p) is given by 
(pi) with (|l5|), where the second-type generating function VF(g, 77) in power series (||) 
is chosen to be (fL6|). 

Remark 3. The convergent radius of the BG-normal form G in power-series (^J) 
vanishes in general (see Moser 1968). In such a case, G is considered only in a formal 
sense. However, like in Remark [l] for K , we will often eliminate 'formal' from such 
formal power series henceforth. 



2.2. The inverse problem 

To define the inverse problem of the BG-normalization appropriately, we review 
the key equation (0) of the ordinary problem from a viewpoint of canonical 
transformations. Let us regard the power series W(q, rf) in (Q) as a third-type 
generating function, a generating function of the 'new' position variables q and 
the 'old' momentum ones 77 (Goldstein 1950), which provides the inverse canonical 
transformation, 

&»7)-(g,p) with £ = -^-p and p = -^Z-p, (18) 
of (§). Equation (@) is rewritten as 

which is combined with ( fL8| ) to show the following. 

Lemma 2.5. Let G(£,7?) of (|) be the BG-normal form for the Hamiltonian K(q,p) of 
(|l|), which satisfies (f|) with a second-type generating function 77) £ imageD gjI) . 
The Hamiltonian K(q, p) is restored from G(£, 77) through the canonical transformation 
( |l8| ) associated with the third-type generating function — W(q, rf) £ imageD q v . 

We are now in a position to pose the inverse problem in the following way: Let the 
Hamiltonian H(q,p) be written in the form, 

n oo 

H (9»P) = 2 E (P* + $ + E ff *(9.P). (2°) 

3=1 fe=3 

where each Hk(q,p) (fc = 3,4, • • ■) is a homogeneous polynomial of degree-fc in (g,p). 
Further, let a third-type generating function S(r], q) be written in the form, 

n oc 
j=l fc=3 



where each S/c(g, 77) (k = 3,4, • • •) is a homogeneous polynomial of degree-fc in (g, 77). 



Definition 2.6 (The inverse problem). For a given BG-normal form, G(£,?7), in power 
series (||), identify all the Hamiltonians H(q,p) in power series (|o|) which satisfy 

H <''-f> = G <^"»- < 22 » 

where the third-type generating function S(q, rj) in power series (^Tj) is chosen to satisfy 
( p2| ) and 

S k (q,r,) G imaged W (fc = 3,4,---). (23) 

We solve the inverse problem in the following way. On equating the homogeneous- 
polynomial part of degree-fc in equation (E2I) is put into the series of equations, 



H k (q, V )-(D q , v S k )(q,T ] )=G k (q,r 1 )-* k (q,r 1 ) (k = 3, 4, • • •), (24) 

where is given by (||). The 9 k (q, rj) is the homogeneous polynomial of degree-fc 
in (q, -q) determined uniquely by H 3 , ■ ■ •, H k _ u G 3 , ■ ■ •, G fe _i, #3) • • given. In 
particular, we have ^3(9,77) = and 



*4(g,7?) 



M 2 V% 



2\d m ) di 3 




2 ^ , (25) 

1 /as 3 V ac, 



Remark 4- As is easily seen from (|9|) and (|25|), ^4 takes a similar form to $4 due to 
(^l|). Indeed, since the substitution W = —S in (Q) provides (g^), for each fc > 4 
can be obtained as & k given by (Q) with H3, • • • , and S3, • • • , S k -i in place of 

K3, • • • , and W3, • • • , W k -i- Such a similarity will be utilized effectively in future 

in writing the program ANFER for the ordinary and the inverse problems in a unified 
form. 



Like in the ordinary problem, we solve (|24j) for H k and S k by using the direct- 
sum decomposition ( [To| ) of V k (q, rj), the vector spaces of homogeneous polynomials 
of degree-fc (k = 3, 4, • • •). Let us decompose H k and ^ k to be 



H k {q,rj) = H^{q,rj) + H^{q,rj), 

^k(q,v) = ^r sc (<i,v) + ^t a '(q,v), 



(26) 



where 



(27) 



K^M^M ekerDW, 
Then on equating kerDq^-part in (p4|), i?£ er is determined to be 

^ ker (g,r ? ) = G fc (g, ?7 )-^ er fe'?). (28) 
Equating imageDg^-part of (|24|), we have 

#r g °(9: »?) - (A^fc) (?> >7) = -*r aE °(9, (29) 
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Since the pair of unidentified polynomials, iJ™ age and Sf. exists in (|29|), 7J™ aEC is not 
determined uniquely in contrast with H^ CI ; such a non- uniqueness is of the very nature 
of the inverse problem. Accordingly, equation ( |29] ) is solved to as 

fffc maE °(<?, v) G imageD^ : chosen arbitrarily, (30) 
S k (q,r)) = (d$> 1 (Hr s %.v)+^r%.v))) M, (31) 
(k = 3, 4, • • •). Now we have the following. 

Theorem 2.7. For a given BG-normal form G(^,ri) in power series the solution 



H(q,p) of the inverse problem is given by (|2C|) subject to (28) and (30), where the 



third-type generating function S(q,rj) in (|22j) is chosen to be ( pl[ ) subject to (31). 
£.3. The degree-25 ordinary and inverse problems 

In the preceding subsections, we have defined the ordinary and the inverse problems of 
BG-normalization, and then find their solutions in power-series. From a practical point 
of view, however, we usually deal with the BG-normal forms not in power-series but in 
polynomial. Indeed, as mentioned in section 1, when we utilize the BG-normalization 
to provide an approximate system for a given system, we truncate the normalized 
Hamiltonian up to a finite degree. Hence it is natural to think of a 'finite-degree 
version' of both the ordinary and the inverse problems (Uwano et al 1999). 

Definition 2.8 (The degree-lb ordinary problem). For a given Hamiltonian K(q,p) of 
the form (|l|) (possibly in polynomial form) and an integer 5 > 2, bring K into the 
polynomial G(£, rf) of degree-25 in BG-normal form which satisfy (Q) up to degree- 
25, where the second-type generating function W{q, n) in (|I) is chosen to be the 
polynomial of degree-25 subject to (||) up to degree-25 and (|7) with k = 3, • • • ,25. 

Remark 5. The reason why we think of only the even-(25-)degree case is that the BG- 
normal form of any 1:1 resonant Hamiltonian consists of even-degree terms only (see 
Remark ||and Lemma 2.2). 



Definition 2.9 (The degree-lb inverse problem). For a given BG-normal form, G(£,?y), 
of degree-25 with an integer 5 > 2, identify all the polynomial-Hamiltonians H(q,p) of 
degree-25 which satisfy (^) up to degree-25, where the third-type generating function 
SJa, ij) is chosen to be the polynomial of degree-25 subject to ( p22| ) up to degree-25 and 
(||) with k = 3, • • • , 25. 

In closing this section, we wish to mention of the way to solve the (degree-25) ordinary 
and the inverse problems on computer algebra. Although the discussion throughout 
this section is mathematically complete, it is not easy to calculate even on computer 
algebra (|l5|,[l6|) and (H|,30,31) as they present because we are faced with a highly 
combinatorial difficulty in calculating $fe and $>k'- To calculate for example, 
G3, • ■ • , Gk-i, S3, ■ ■ • , Sk-i, and H3, • ■ • , Hk-i have to be kept on computer. What 
is worse is that the dimension, J2h=i Od-i)' °^ to which Gg, Sg and Hg belong 
rises in a combinatorial manner as I increases, where the symbol Q indicates the 
binomial coefficient. These facts will cause a combinatorial increase of the memory- 
size on computer required for calculation. To get rid of such a difficulty, we break the 
transformations, (|4|) and ( p2| ) into a recursion of certain canonical transformations of 
simpler form, which will be presented in Appendix for the degree-4 case. 
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3. Example: The one-parameter Henon-Heiles system 

In this section, we take the one-parameter Henon-Heiles-Hamiltonian 

2 

2 



1 2 

^.(9>P) = 5£(P? + 9?)+9i92+W2 (MGR.), (32) 

as an example to illustrate how the ordinary and the inverse problems are proceeded. 
As is mentioned in section 1, this example inspired the present work. 

3.1. The degree-A ordinary and inverse problems 

In order to present the results in compact forms, the complex variables defined by 

Zj = qj + ipj, Cj = £j + iVj 1? = 1 > 2 )> ( 33 ) 
will be of great use. In terms of z and z, is written in the form, 

K^(q,p) = -(ziZi + z 2 z 2 ) + ^(z 2 + 3z^z 2 + 3z 2 zj + z\) 

2 x 8 (34) 

+ -{z\z 2 + zfz 2 + ~Z\z\ + ~z\~z\ + 2z{ziz 2 + 2z{z{z 2 ). 
8 

Using a prototype of ANFER (Uwano et al 1999, Uwano 2000), we see that the BG- 
normal form for is given, up to degree-4, by, 

G^,J7) = i(CiCi+C2C 2 ) 

+ ^ {-5Ci 2 Ci - 45 M 2 C|C2 - (8 + 36/.)CiC2CiC 2 ( 35 ) 
+3 M Ci 2 C2 + ZKlt - 6Ci 2 ?2 - 6C|Ci} , 

which is well known (see Kummer 1976, Cushman 1982). We proceed to the 
degree-4 inverse problem for the BG-normal form G M in turn. Solving ( p2| ) with 
in place of G by ANFER, we have the Hamiltonians in polynomial of degree-4 of the 
following form as the solution; 

2 

H^(q,p) = \ J2(p 2 j + qf) + H^ 3 (q,p) + JT„, 4 (<?,p) (36) 



2 



with 



H^,3(q,p) 



aizf + a 2 z\z 2 + a 3 zxz 2 + a^z\ + a^z\z\ 

„2- , „ „ „ — , „ „ „ -= | „ 2— , „ _2- 



a 6 z 1 z 2 + a-?ziz 2 zi + a s ziz 2 z 2 + a 9 z 2 zi + a w z 2 z 2 (37) 



+ a\z\ + a 2 z\z 2 + azz\z\ + a±z 2 + a$z\z% 

+ a & z 2 z\ + ~a~rZiZiz 2 + a s z 2 z{z 2 + agz{z 2 + ai z 2 z 2 



and 



C\z\ + C 2 z\z 2 + C-sz\z\ + CiZiZ 2 + c 5 z 2 + c 6 zfzi + cjz\z 2 
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+ c$zjz 2 z 1 + c 9 z\z 2 z 2 + c X qZ\Z%Z\ + c u ziz%z 2 + c^zf^i + c 13 z 2 z 2 
+ c\~z\ + c 2 ~z z {z 2 + c^z 2 {z% + C4Z1Z2 + c 5 z 2 + c 6 zizl + c-jz 2 ~z\ 

+ C$Zl~zl~Z 2 + CgZ 2 zlz 2 + C W Z{ZlZ% + Cn222l2 2 + C\ 2 Z{Z 2 + Cl 3 Z 2 Z?, 

+ 8 (09010^2^1^2 + a 9 a 9 ziz 2 ziz 2 + a\oa 9 ziz 2 z% 
+a 5 a e ziz 2 zl + a 6 a 5 zfziz 2 + a 6 a 6 ziz 2 Ziz 2 ) 
a\5,\z{z x + a\ 9 ~a\ 9 z 2 ~z 2 + a, x a^z 2 'z 2 + a^a^z^z^ 

+ 4 (a s a 5 ziz 2 ziz 2 + a$a e z 2 ziz 2 + a$a 9 z 2 zl + a 9 a 7 z\z 2 ~z\ 

+ agasz 2 ~Zi + a\H 2 z 2 z~\~z 2 + aioa 7 ziz 2 zi~z 2 + a 2 aiziz 2 z 2 
+ 03042122^2 + 0403^2^122 + a 5 a 8 ziZ2ZiZ2 + 06072^2 (38) 
+a6a82i2;2^2 + 07010^1^2^122 + ajaez 2 ! 2 + 07092^2122) 
8 , _ _ _ 

+ g (0 2 a 2 2l222l22 + a 3 a 3 2iZ22l22) 

+ 2 (a 8 aioZiZ222 + a8072i2i2 2 + a 8 a 7 2i2222 + a s a s z 2 zl 
+ a\a 3 z\z\ + a\Qa%z^z\z 2 + a 2 a4Z 2 ~z 2 + a 3 a\z 2 'z\ 
+ ai<i 2 z 2 ~z\ + 05072^2122 + 070521222^ + 07072^2^ 

+ a 7 a s ziz 2 zl + a 7 a s z 2 ziz 2 - a%a§z\z 2 ~z\ - ara 9 ziz 2 zl 

2 2 2 2 2 2 2 

— aga^,z 2 z 1 — a 9 a 7 z 2 z\z 2 — a 9 a 9 z 2 z 2 — a\ 9 a§z 2 z x 

2 2 2 2 2 2 2 \ 

—a^a 9 z 1 z 2 — aeaioz 1 z 2 — aeaez 1 z 1 — a§a%z 1 z\z 2 J 

+ g (+02032^2122 + 2 a 3 2l2222 + 3 a 2 2l22Zi + a 3 a 2 2f Zl2 2 ) 

+ I (a 2 a 2 z 2 z\ + a 3 a 3 z 2 zl) 

+ — (-8z 1 z 2 z 1 'z 2 ~ hz\~z\ - &z 2 ~z\ — 6z 2 z 2 
48 

— 36/i2iZ22i22 — 45n 2 z 2 z 2 + 3/iz 2 z 2 + 3/IZ 2 Zi^ , 

where (h = 1, ■ • • , 10) and q (£ = 1, • ■ • , 13) are complex- valued parameters chosen 
arbitrarily. In (38), the polynomial with the coefficients (ct) expresses H 1 ™^ and the 
polynomial whose coefficients are written in terms of (ah) and (/&) does H^f\. From 



the lengthy expression (36, 37, 38), one might understand an effectiveness of computer 
algebra in the inverse problem. 

It is worth pointing out that if we choose (ah) and (ci) to be 

ai = a 3 = «5 = «8 = a 9 — 0, 

1 3 M (39) 

2a 2 = 2a 6 = a 7 = -, 3a 4 = aio = — , 

and 

Q = (£= 1, ■ • • , 13), (40) 
respectively, becomes equal to the one-parameter Henon-Heiles Hamiltonian K^. 

3.2. The Bertrand-Darboux integrability condition 

We wish to find the condition for (ah) and (ci) to bring into the Hamiltonian of 
the perturbed harmonic oscillator with a homogeneous-quartic polynomial potential 
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(PHOQP). To bring H„ into a PHOQP-Hamiltonian, we have to make -ff M ,3 vanish. 
The vanishment of is realized by the substitution 

a h = (/i = l,---,10) (41) 

in (37). We bring i? M ,4 1 a =o (<7, p) into a homogeneous polynomial of degree-4 in q in 
turn. To do this, we have to find the set of non- vanishing ((q), (A m )) for which the 
following identities for z hold true; 

Ai<?i =ciz t + c 6 zfz 1 - —zfzl + c{z[ + c§z{z%, 
4o 

A2gi?2 = c 2 z\z 2 + +c 7 zfz 2 + c 8 zlz 2 zi (42) 

+ C 2 zfz 2 + C-jZ 2 ~z\ + CgZiZ{Z2, 

^qiq 2 = C11Z1Z2Z2 + ci 2 Z2~zi + cazizI 

+ C4J1ZI +ci 2 z{z% + cnz 2 zizl, 
15 

A 4 <72 = c $ z 2 + azzl~z 2 - —[i 2 z\z\ + ci 3 z 2 zl + c 5 z^, (43) 

2 2 2 2 2 2 

A5<7l<72 = C\qZ\Z^Z\ + c 32l^2 + C9Z1Z2Z2 

+ c w z{z{z1 + c 3 zfzl + c 9 z 2 z~lz 2 

/-i 2 2 M 2 2 

H Zi Zn H ZnZ, z-\z 2 z\z 2 

16 1 2 16 2 1 4 

1 1 9 2 1 2 2 

- -^z x z 2 ziz 2 - -^z 2 - g^ 2 i' 

where A m (m = 1, • • • , 5) are real- valued parameters, and qj — (zj + z~j)/2 (j = 1,2). 
To make the identities hold true from the first to the fourth, we have to choose the 
parameters appearing in those identities to be 

5 

Ai = 16ci = 4c 6 = - — , A 2 = c 2 = c 7 = c 8 = 0, 

f>? (44) 
A 3 = c 4 = cu = C12 =0, A 4 = 16c 5 = 4ci 3 = — — . 

The fifth identity holds true if and only if the overdetermined system of equations, 
As = 16c 3 = 8c 9 = 8cio, 

A 5 = -3 M -| (45) 
A 5 = M - 2, 

admits a solution. By a simple calculation, we see that ( f45| ) admits the solution, 

A 5 = 16c 3 = 8c 9 = 8c 10 = -~, (46) 
if and only if fi satisfies 

*-\. (47) 
To summarize, we have the following. 

Theorem 3.1. The one-parameter Henon-Heiles Hamiltonian K^(q,p) shares its BG- 
normal form up to degree-4 with the Hamiltonian of the perturbed harmonic oscillator 
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with a homogeneous-quartic polynomial potential (PHOQP) if and only if the 
parameter fj, satisfies (fl7|). Under (^), the PHOQP-Hamiltonian sharing the BG- 
normal form with -K1/3 is given by 

Qb>p) = \ YM + - + 6 <^f + (48) 



After theorem 3.1, one might ask whether K 1 / 3 and Q have specific meanings or not. 
Those who are familiar with the separability of dynamical systems may recognize 
immediately that -ffi/3 and Q are separable in q\ ± q 2 (see Perelomov 1990). We 
can hence answer this question affirmatively owing to the Bertrand-Darboux theorem 
concerning not only the separability but also the integrability. The theorem is stated as 
follows (see Marshall and Wojciechowski 1988, Yamaguchi and Nambu 1998, Groschc 
ct al 1995, for example). 



Theorem 3.2 (Bertrand-Darboux) . Let F be a natural Hamiltonian of the form, 

2 

2 



1 2 

F(q,p) = -Y,P 2 3 +V(q), (49) 
j=l 

where V(q) a differentiable function in q. Then, the following three statements are 
equivalent for the Hamiltonian system with F. 

(1) There exists a set of non- vanishing real-valued constants, (a, j3, j3' , 7, 7'), for which 
V(q) satisfies 



fd 2 V d 2 V\ , 

\~dqj ~ ~dqj ) (~ 2aqiq2 -I 3 l2~f3qi+ 7) 



i>2] ' (aql-jqf + fa-P'qi+l') (50) 



dq 1 dq 2 

dV dV 
+ T—(6aq2 + 3/9) - t— (6agi + 3/3') = 0. 
oqi dq 2 

(2) The Hamiltonian system with F admits an integral of motion quadratic in 
momenta. 

(3) The Hamiltonian F is separable in either Cartesian, polar, parabolic or elliptic 
coordinates. 



Due to the statement (2) in theorem |3.2| , a natural Hamiltonian system with F is 
always integrable if ( ]50| ) holds true. In this regard, we will refer to ( |50| ) as the 
Bertrand-Darboux integrability condition (BDIC) henceforth. 

For the PHOCP's and the PHOQP's, Yamaguchi and Nambu (1998) have given a 
more explicit expression of the BDIC ( |50| ) convenient for our purpose: 

Lemma 3.3. Let J-^ k \q,p) (k = 3,4) be the Hamiltonians of the form, 

F {k) [q,p) = \jZ(p] +<$) +^ k) {q) (k = 3,4), (51) 
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with 



V i3) (q) = hq\ + hq\q2 + hqidl + Uql (52) 



P (4) (?) = giqf + 9iq\qi + a?,q\q\ + gmq% + g$qt, (53) 

where fk (k = 1, ■••,4) and gg (I = l,---,5) are real-valued parameters. For 
the PHOCP with P^(q,p), the BDIC © is equivalent to either of the following 
conditions, (54a), (546), or (54c); 

3(/i/ 3 + hU) - {fl + fl) = 0, (54a) 

ft = 2/ 3 , h = h = 0, (546) 

h = 2/ 2 , fi = h = 0. (54c) 



For the PHOQP with T^{q,p), the BDIC © is equivalent to either of the following 
conditions, (55a) or (556); 



(55 a) 



.93 = 2#i = 255, 
92 = 9i = 0, 

9ff! + 4. 9 | - 24 5lff3 - 9 3 234 = 0, 

9.9! + 4<?§ - 24 53 55 ~ 9 5 254 = 0, (556) 
(92 + 94)93 ~ 6(31.94 + 5235) = 0. 

The integrability of the Henon-Heiles system with /i = 1/3 and the PHOQP with Q 
can be confirmed now by showing that ( |50| ) holds true both for Ki/ 3 and Q. Let us 
start with showing the integrability of the one-parameter Henon-Heiles system with 
li = 1/3. If we choose (fh) to be 

/i = / 3 = 0, f 2 = l, / 4 = i (56) 

the Hamiltonian J 7 ^ becomes ifi/3. Evidently, the (//J given by ( |56| ) satisfy 
the BDIC (54a), so that the one-parameter Henon-Heiles system with \x = 1/3 is 
integrable. 

We proceed to showing the integrability of the PHOQP with Q given by ( ^8|) in 
turn. In order to bring J 7 ^ into Q, we choose (gg) to be 

5 

6.9i = 93 = Q95 = 92 = 9i = 0. (57) 

By a simple calculation, we see that the (gg) given by ( |57|) satisfy the BDIC (556), so 
that the PHOQP with Q is integrable. To summarize, we have the following. 

Theorem 3.4- If the one-parameter Henon-Heiles system and the perturbed harmonic 
oscillator with a homogeneous-quartic polynomial potential admit the same BG- 
normalization up to degree-4, then both dynamical systems are integrable in the sense 
that they satisfy the Bertrand-Darboux integrability condition. 



4. Extension: The perturbed harmonic oscillators with 
homogeneous-cubic polynomial potentials 
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The degree-A ordinary and the inverse problems of PHOCP 

In this section, the discussion in section 3 for the one-parameter Henon-Heiles system 
is extended to the perturbed harmonic oscillators with homogeneous-cubic polynomial 
potentials (PHOCP 's). To start with, we solve the degree-4 ordinary problem for the 
PHOCP-Hamiltonian T^\q,p) given by @ with (||). By ANFER, we see that the 
BG-normal form for JF( 3 ) is given, up to degree-4, by 

= 5 (C1C1 + C2C2) 

~ §(/iC?c5 + flClcl) - -(A/sCiCaCiCa + /2/4C1C2C1C2) 

- |(/i/ 2 Ci 2 CiC 2 + /1/2C1C2C1 + /3/4C1C2C2 + /3.AC2C1C2) 

- ^(/ 2 /3Ci 2 CiC 2 + /2M1C2C1 + /2/3C1C2C2 + /2/ 3 C 2 2 CiC 2 ) (58) 

- ~(/2&C»CiC2 + /IC1C2C1C2) - J(/fCi 2 Ci + /!C 2 2 C2) 

- lifklcl + ficicl + mil + ficlcl) 

+ ^(/l/3Cl 2 C2 + /1M2C1 + /2.AC1C2 + hMltl). 

Note that if we choose (fh) to be 

/i = / 3 = 0, / 2 = 1, h = n, (59) 

Q becomes G^, the BG-normal form for the one-parameter Henon-Heiles Hamiltonian. 

We solve the degree-4 inverse problem for the BG-normal form Q in turn: By 
ANFER, we have the following polynomial of degree-4 as the solution; 

1 2 

U ^P) =nE (Pi + if) + U ^P) + (60) 



with 



2 



7~(-3(q,p) = ai^i + a 2 z\z 2 + a-zz\z\ + a^z\ + a 5 zfzi + a 6 ZiZ 2 
+ a 7 ziz 2 zi + a s ziz 2 z 2 + agz^zi + a w zlz 2 
+ ai~z\ + a 2 z\z 2 + az~zi~z 2 + a±z\ + a^z{z\ + a 6 z 2 z~l 
+ a 7 ziziz 2 + a 8 z 2 ziz 2 + a<jZ-{z\ + a w z 2 ~z 2 , 



(61) 



and 



Hi(q,p) 



= c\z y + c 2 z\z 2 + czz\z\ + c 4 ziz| 
+ C5Z2 + c e zl~zi + c 7 z\~z 2 + c$zlz 2 Zi 
+ co,z\z 2 ~z 2 + C10Z1Z2Z1 -I- c\\z\z?£z 2 + ci-izfzi + ci 3 z|z 2 
+ c{z\ + ~c 2 ~z\z 2 + ~c^z\z\ + CiZ\z% 
+ C5Z2 + c 6 zizl + c 7 z 2 zl + c 9 ,zi~z\~z 2 

+ CgZ 2 zlz 2 + CloZlZlzl + Cll 22^1^2 + C1221Z2 + Cl3Z 2 Z?, 

+ 8(aea6Ziz 2 ziz 2 + agagziz 2 ziz 2 + a<^a§z\z 2 ~z\ 
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, — 2 , — 2 , — — 2\ 

t aQa 5 z 1 ziz 2 + agai z 2 ZiZ2 + aioagZxZ2Z 2 ) 

/ \ 

+ &{a\a\z 1 z 1 + a±a^z 2 z 2 + a^a^z 1 z 1 + a\Qa\QZ 2 z 2 ) 

+ 4:(aia 2 zfziz 2 + a%agz{z 2 + a 3 a i ziz 2 zl + a 5 a s ziz 2 ziz 2 

+ a§aiz{z 2 + a e a 8 ziz 2 zl + aiagz\~z-{Z2 + arai ziz 2 ziz 2 ) 
+ 4:(a,2aiziz 2 zl + aia^z^xzi + a%a$z\z 2 ~z\z 2 + aiaez 2 zl 

+ asaez 2 ziZ2 + ag~a-!Z\Z2~z\ + awajzi Z2~z\zi + aga&ZyZx) 

8, _ _ _ 

+ -^\a2a2Ziz 2 ziz 2 + a 3 a 3 ziz 2 z 1 z 2 ) 

i <">/ — 2—2 i — 2—2 i — 2—2 — 2— 2\ 

+ 2(— a§a§z l z 1 + a-;a'jz 1 z 1 + a 8 a 8 z 2 z 2 — agagz 2 z 2 ) (o2) 

. 2 o «= 2 2 — 2 — — 2 2 

+ 2(a\a 3 z 1 z 2 + a20AZ 1 z 2 + a^aiz 1 z\Z2 — a$agz 1 z 2 

— a§a%z\z\z 2 — a§a~\gz\z\ + aia^,z\z 2 z\ 

+ aia%z^z\Z2 — a7~agZiZ2~z\ + a 8 a\§z\z 2 £ 2 ) 
+ 2(a 3 aiz 2 z 1 + a4& 2 z 2 ~z 1 + a-ia^z\Z2Z\ — aga$z 2 'z 1 

- a$a§ziZ2~z\ — a 10 a 6 z%'z~l + a^a^z^z^ 

+ a 8 a7ZiZ2zl - a 9 a7Z 2 ziz 2 + a w a 8 z 2 ziz 2 ) 

+ ^(0203^1^1^2 + a 2 a 3 z 1 Z2'z~l + a 3 a 2 ziz 2 z\ + a^z^xz^ 
+ \ {p,2a2z\'z\ + a 3 a 3 z 2 zl) 

~ t|(/i Z 1 7 1 + fi Z 2^l) - ^(flf3ZlZ2ZlZ 2 + hf^\ziz\z%) 
lb 4 



- -g{flhzlziZ2 + flf2ZlZlz\ + f 3 fiZlZ 2 zl + j 3 f4,zlziZ2) 

'-{zfzxZ2 + ziz 2 zl + ziz 2 zl + z 2 ziz 2 ) 
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- - (fiz 1 z 2 z 1 z 2 f 3 z 1 z 2 z 1 z 2 ) 

~ g(/2 z l z 2 + /l z 2 z l + .fi z l z 2 + fi z 2 z l) ~ ^(/2 Z 1 Z 1 + fi z 2 z 2) 

+ ^{fihztz\ + fif 3 zM + hhz\z\ + hf^\z\), 

where ah (h = 1, • • • , 10) and ci (£ — 1, • • • , 13) are the complex-valued parameters 
chosen arbitrarily, and (k = 1, •••,4) the real-valued parameters in V^\q) (see 
([52|)). Namely, we have 46-degree-of freedom in the solution, TL, of the inverse problem 
of the PHOCP with (/&) fixed. Note that if (fk) are chosen to be ( |59|) then TL becomes 
H u given by (^6|)-(38). Further, if (ah), (ci) and (fu) are chosen to be (39), (|40| ) and 
([59]), respectively, then TL becomes given by (p2|). After (|60|)-(62), one might 
understand more than after (|36|)-(38) the necessity of computer algebra in the inverse 
problem. 



4-- 2. The Bertrand-Darboux integrability condition 

We wish to find the condition for (ah), (q) and (fk) to bring TL into the PHOQP- 
Hamiltonian J 7 ^ (q, p) defined by ( pl|) w ith j5^ ) . In order to make TL 3 (q, p) vanish (see 
(|6l|)), we have to choose (ah) to be (flip. Our next task is then to bring TL4,\ a =o(Q,p) 
into a certain V^(q). Let v n (n = 1, ■ • • , 5) be real- valued parameters. Then we see 



C\2Z\Z% + C\\ZlZ\z\ + CAZ\Z% 
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that 7-^4 | a= o takes the form V^\q) if and only if the following identities of z hold true 
for a non-vanishing ((ch), iy n )); 

A A 3 45 f? ~t~ 5 /o 9 2 — 3 — 4 

v\q x = c\z x + c 6 z\zi — z y z x + c%Z\Z x + az x , 

48 

4 4 , 3- 45/| + 5/| 2 _ 2 _ 3 4 

^2<? 2 = C 5Z 2 + c 13 z 2 z 2 — z 2 z 2 + C 13 Z 2 Z 2 + C 5 Z X , 

V3qf <?2 = c 2 z\z 2 + c 7 zfz 2 + cszfz^x 

+ ~C2~z\z~2 + C7zfz2 + CsZ 2 {Z2Z\ 

5 

- ^( 3 /l/2 + hh){zlziZ 2 + ZlZ 2 z\), 

^4(7192 = C4Z X Z 2 + C 1 lZ 1 Z 2 Z 2 + C\2Z 2 Zl (63) 

C12 
5 

- ^( 3 /3/4 + f2h){zlziZ2 + ZlZ2Z%), 

v$q[q\ = c z z\z\ + c 9 zfz 2 z2 + c w zizl~z\ 

+ C10Z1Z1Z2 + CgZ2zlz2 + C3zfzl 

+ 3^/1/3 + hh - 2/| - 2/fX^ +z\zl) 

~ ^(9/i/3 + 9/2/4 + 2/ 2 2 + 2/ 2 )z 1 Z 2 ^2, 

where v n (n — 1, • ■ • , 5) are real- valued parameters and Qj = {zj +~Ej)/2 (J = 1, 2). 
As a necessary and sufficient condition for the identities to hold true from the first to 
the fourth, we have 

^(9A 2 + /2 2 ), 

yC 8 = -^/ 2 (3/ 1 + / 3 ), (64) 

yen = -y/ 3 (3/ 4 + h)- 

As for the fifth identity, we see that it holds true if and only if the following 
overdetermined system of equations, 

^5 = 16c 3 = 8c 9 = 8ci , 

^5 = (/i/3 + /2/4)- 2(/ 2 2 + / 3 2 ), (65) 

^5 = -3(/l/3 + /2/4)-^(/ 2 2 + /3 2 ), 

admits a solution. Surprisingly, the Bertrand-Darboux integrability condition (54a) 
for PHOCP's comes out as a necessary and sufficient condition for ( |65| ) to admit a 
solution ! Under (54a), equation (65) admits the solution, 

v b = 16c 3 - 8c 9 = 8c 10 = (/1/3 + /2/1) - 2(/| + fl), (66) 



V\ 


= 16ci 


= 4c 6 = - 


l>2 


= 16c 5 


= 4c i3 = - 


t>3 


= 16c 2 


= 16c 7 = 


Vi 


= I6C4 


= I6C12 = 
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which is combined with ([ll]) and (64) to bring TL into the PHOQP-Hamiltonian, 

- ^OA 2 + fM - y(3/i + f 3 )hqU2 ~ l(fl + fM<& (6?) 

10 ^ 

- y (3/ 4 + / 2 )W - ^(9/1 + / 3 2 )<Z 2 4 , 

where (fh) are subject to the BDIC (54a). We note here that the coefficient of qfq% in 
( |67| ) is obtained by combining ([36]) with (54a). To summarize, we have the following. 

Theorem J^.l. The perturbed harmonic-oscillator Hamiltonian J-^ with a homoge- 
neous cubic polynomial potential shares its BG-normal form with the perturbed 
harmonic-oscillator Hamiltonian J-^ with a homogeneous-quartic potential up to 
degree-4 if and only if the PHOCP-Hamiltonian satisfies the Bertrand-Darboux 
integrability condition (54a). Under (54a), the PHOQP-Hamiltonian J 7 ^ sharing its 
BG-normal form with that PHOCP-Hamiltonian J 7 ^ is equal to Q given by (|67j). 

We are now in a position to show the integrability of the PHOQP with Q subject to 
(54a). It is easy to see that Q is given by ( [5l| ) with j53| ) under the substitution, 

5 10 5 

.91 = -^(9/? + /!), <72 = - ¥ (3/i + / 3 )/2, 9 3 = -^fi+m, 

1 5 

<?4 = -y(3/4 + / 2 )/3, 55 = -^(9/| + / 3 2 ). (68) 

A long but straightforward calculation shows that (gi) given by (68) with (54a) satisfy 
the BDIC (556), so that the PHOQP with Q is integrable. 

Theorem 4-2. If the perturbed harmonic oscillators with a homogeneous-cubic 
polynomial potential and with a homogeneous-quartic polynomial potential share the 
same BG-normal form up to degree-4, then both oscillators are integrable in the sense 
that they satisfy the Bertrand-Darboux integrability condition. 



5. Concluding Remarks 



As is shown in the previous section, the new deep relati on b etwe en the BDIC for 
PHOCP's and that for PHOQP's are found (see theorems |t| and ||J). The results 
obtained in the present paper are expected to provide several interesting subjects, 
which are listed below. 



(1) A further generalization of theorems 4.1 and 4.2 will be worth studying: A 



conjecture is posed as follows, which is investigated now. 

Conjecture 5.1. Let (q,p) and J^i 2r - 2 ) (g,p) (r = 

perturbed harmonic oscillator Hamiltonians of the form 

^ ( %,P) = ~X>i + $ + ^ (r) (<7), 
F (2r - 2) (q,p) = \J2(P" +® + V (2r - 2) (q), 



3,4,- 



denote the 

(69) 
(70) 
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where V^ 2r 2 ) is a homogeneous-polynomial potential of degree- (2r — 2) and 
of degree-r subject to 

V {r] {q) e imageDW. (71) 

If the perturbed harmonic oscillators with the Hamiltonian J 7 ^ and with J^i 2r ^ 2 ) 
share a BGNF up to degree- (2r — 2) then both oscillators are integrable in the 
sense that they satisfy the Bertrand-Darboux integrability condition. 



(2) As is pointed out after theorem 3.1, the Hamiltonians, K\i^ and Q have a 
significant feature other than integrability; they are separable in q± ± q%. Such 
a separability can be found also in !F^ 3 ' with (54a) and Q, the generalization of 



K 1 /z and Q, respectively. In fact, we can find, from theorem 4.1, the PHOCP's 
with 

P< 3 >(?) =a( qi +q 2 ) 3 + b( qi - q 2 f (a, b € R), (72) 

as special cases of PHOCP's subject to the BDIC (54a), which cover all the 
PHOCP-Hamiltonians separable in qi ± q 2 (see Perelomov 1990). Further, 



theorem 4.1 implies that each of the separable PHOCP's (subject to (f72j)) shares 



the same BG-normal form with the PHOQP with 

V^iq) = -5{(a a (<fc + q 2 f + b 2 ( Ql - a 2 ) 4 }, (73) 

which is also known to be separable in qi ± q 2 (Perelomov 1990). It is worth 
noting that all the PHOQP 's with (j73) cover a quarter of the class of all the 
PHOQP-Hamiltonians separable in q\ ± q 2 in view of a 2 , b 2 > 0. Thus we reach 



through theorem 4.1 to the four-to-one correspondence between the PHOCP's 
with (f?|) and PHOQP 's with @ separable in qi ± q 2 . Since with (54a) 
and Q are thought to include several classes of Hamiltonians separable in several 
coordinate systems other than q\ ± q 2 , the separability will be worth studying 
extensively from the BG-normalization viewpoint in future. 



(3) The perturbed oscillators referred to in theorems 4.1 and 4.2 are expected 
to provide good examples of the quantum bifurcation in the BG-normalized 
Hamiltonian systems (Uwano 1994, 1995, 1998, 1999): Since the perturbed 



oscillators referred to in theorems 4.1 and |4.2j are integrable, their quantum 



spectra are expected to be obtained exactly. Then we will be able to think 
of whether or not the quantum bifurcation in the BG-normalized Hamiltonian 
system for those oscillators approximates the bifurcation in these oscillators in a 
good extent. 

(4) In section 1, two approaches to the BDIC (or BD-theorem) have been mentioned. 
As another approach to the BDIC, the work of Yamaguchi and Nambu (1998) 
is worth pointing out, in which the BDIC for the PHOCP's and the PHOQP's 



(see Lemma 3.3) came out from the renormalization of Hamiltonian equations. 
It will be an interesting problem in future to study a relation between the BG- 
normalization and the renormalization. 
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(5) As is easily seen, the solution @ |l[ 62) of the inverse problem for Q admits fifty 
real-valued parameters. We may hence expect to obtain other integrable systems, 
so-called the electromagnetic type (Hietarinta 1987), for example. 

On closing this section, we wish to mention of the role of computer algebra in the 
present paper: Without computer algebra, it would have been very difficult to find 
theorems 3.1, B.4, 4.1 and \ 
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Appendix A. Mathematical basis for coding ANFER 

In this appendix, we describe in mathematical terminology the algorithm of ANFER, 
a symbolic computing programme, for the degree-4 inverse problem on REDUCE 3.6 
coded by the author. Although only a primitive prototype exists, those who 
are interested in ANFER will be able to see its source-code at Uwano (2000), 



http://yang. amp. i. kyoto-u.ac.jp/~uwano/. After the algorithm, we present a key 



lemma supporting the algorithm. The general-degree case has been reported in Uwano 
et al (1999) briefly, and will be discussed in detail in a future paper by the author. 
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Appendix A.l. Setting-up 

The setting-up of ANFER is made as follows for the degree-25 inverse problem. Let 
SCO^CO, (h = 3 7 . . . , 28) be the third-type generating functions of the form 

SW(f W jjf- 1 )) = - ^ ^ty^-S^W r^" 1 *) (ft = 3, • • • , 25),(Al) 

where each 5ft are the homogeneous polynomial of degree-ft. With 5^, we associate 
the canonical transformations, 



ith 



ru : {^ h - l \r 1 {h - 1) ) - (f^W ) (ft = 3, ■ • .,25), 



(h-i) 



drj (h- 



T) 



(ft = 3, •••,25). 



(A2) 



(A3) 



From a given BG-normal form G of the form (ph, we define the initial Hamiltonian 
if ( 2 ) to be the degree-25 polynomial form, 

2 25 



j'=i 



fc=3 



In ANFER, the degree-25 inverse problem, 

H {h) ( Ah) dsih) \ _ ff (h-i) ( dsW „(h-i) } 
SW(£W „W) e unageDgw,^-!) with @, 



(A5) 
(A6) 



for ifC 1 " 1 ) are solved recursively with ft = 3, • • • , 25. The resultant ij( 25 ) provides the 
solution H of the inverse problem for the BG-normal form G (Uwano et al 1999). 

Appendix A. 2. Algorithm for the degree-A inverse problem 

We show how (A5, A6) is solved in the degree-4 case: 
Step-I: Solving (A5, A6) with ft = 3 

Equating the homogeneous-parts of degree-2, 3 and 4 in (A5) with ft = 3, we have 



fff(£ (3 W 2) ) = ffi 2) (£ (3 W 2) ), 

flf tt (3) , * (9) ) - ,„«> 5 3 ) (c (3) ,^ (2) ) = # f (c (3) , ^ ) 



fff(e (3) ^ (2) )-E 



1 / dS; 



dS 3 dH. 



(3) 



(3) 



(A7) 
(A8) 

(A9) 



(J(3) )?J (2)) 



1 dSa_ 
5 Uf, 



95, flff. 



(2) 



(t (3) ,v (2) )/ 
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In a similar way to (p4|-31), (A8) is solved to be 

#(3)kcr (£(3)^(2)) = ^(2)^^(3)^(2)) = ^ 

(i?^ 3) ) imagc (f (3) ,77 (2) ) e imageD^ (2) : chosen arbitrarily, 
53^(3)^(2)) ' ' (A10) 



£)(3) 
?,») 



-1 



The H[ 3) is given by (A9) with the substitution (|A10|) into i^ 3) and <S 3 . 



Step-II: Solving (A5, A6) with h = A 

Equating the homogeneous-parts of degree-2, 3 and 4 in (A5) with h = 4, we have 

^^(4)^(3)^^(3)^(4)^(3)) (fc = 2j3)j (AH) 

^(C( 4 ),^)-(^, i ,( 3) 5 4 )(^,7 ? ( 3 ))=ili 3 \^),7 ? ( 3 )). (A12) 
In a similar way to (p4[-31), (A12) is solved to be 

ff ( 4 )ker^(4) )J7 (3)) =jff (4)ker (e (4) )?? (3)) 

i/| 4)imago (^ (4) ,77 (3) ) € imageD^i, (3) : chosen arbitrarily, 

5 4 «W,77 (3) ) ' (A13) 

= (^[^, K^k*,,) -tfi 3W l( 9 ,n))) « (4) ^ (3) )- 

We are now in a position to show that H^{q,p) and S(q,n) = — Y2j=iQj r lj ~ 
Sfe=3 v) are identical with H(q,p) and S(q,T]), respectively, the solution of the 
degree-4 inverse problem for G. Equations (|25| ), Gz{q,n) — and ^ 3 (q,rj) — (see 
the line above (|25|)) are put together with (|28|, 30, 31) to show that Hs(q,p) is equal 
to (q,p), so that we have Ss(q,T]) = S3(q,n). Under H 3 (q,p) = (q,p) and 
Sz{q^ v) — ^(q, f]), we see that the second term in the r.h.s. of (A9) is equal to — \&4 
given by (|5j) with (q,ri) — (^ 3 \?/ 2 * ) ). This fact and (A12) are put together with 
( |Al'j| ) to show that H^ 4 \q,p) and S±(q,r]) are identical with H±(q,p) and S^q^rj), 
respectively. 

What is characteristic of the algorithm above is that the procedure of solving ( pi]) 
with k = 4 is divided into a pair of steps, Step-I and Step-II (extending to the degree- 
25 case, we have 26 — 2 steps). Although it is not so significant in the inverse problem 
of lower-degree (like degree-4), such a division will contribute a lot to the reduction 
of the memory-size required through computation. 

Appendix A. 3. Composition of canonical transformations 

We wish to give an explicit expression of the generating function associated with 
the composition T4 o T3 of canonical transformations, T3 and T4, which support 
mathematically the algorithm given above. The general-degree version can be found 
in Uwano et al (1999) without a proof. Such an expression is really important because 
little is known explicitly of the composition of non-infinitesimal canonical although 
well-known is that of infinitesimal ones (Goldstein 1950). 
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Lemma A.l The composition T4 o t 3 : (£( 2 ) , rfi*' ) — > {^^iV ) °f the canonical 
transformations T3 and T4 defined by ( |A2[ ) is associated with the third-type generating 
function 

S(£ (4 W 2) ) -E^+ 5(3) (C^ (2) )+5 (4) (e (4) ^), (A14) 
where f and 77 are the functions of and r]^ uniquely determined to satisfy 



around (£ (4) ,?7 (2) ) = 0. The 5 admits the expansion, 

2 

<5(g,»7) = - ^QjVj - s 3(q,v) -S 4 (q,v) +04(9,77), 
J'=l 



04(9,77) 



(, 



v £(«?+tf)-°). 

\ i=i 



Proof: Since the third-type generating functions, iS'^'s, satisfy 
_ € (/.-i) dtJ (fc-i) _ qWdgW = dsW (h = 3, 4) 

(Goldstein 1950), we have 

-e (2) d/? (2) - 77 (4) d£ (4) 

= 4 (£ $° •tf ) + 5(3) ^ (3) < ^ (2) ) + 5(4) « (4) ' ^ (3) ) 



(A15) 

(A16) 
(A17) 

(A18) 

(A19) 



We express £( 3 ) and 77^' in terms of £ and 77 in turn: Applying the implicit function 
theorem (Spivak 1965) to the map, 



(£ (3 W 3 W 2 U (4) )eR 



(A20) 



around (£^, 77^, 77^, £( 4 )) = 0, we find the unique pair of functions, £ and 77, in £( 4 ) 
and 77(2' which satisfy t(£, 77, 77^), £<4>) = around (£ (4) ,t/ 2 )) = 0. This show s (|Al^ ) 
with d A 1 5|) . We move on to the proof of (A16) with (A17) in turn. Equation ( A15[) is 
expanded in (£ < - 4 - ) , 7/ 2 ') as 



8Sa 



where o 3 and o 3 ' indicate the terms of the order higher than 3 in \J^2'j = i(£,j' ±>2 + Vj 
Substituting (A21) into (A14), we have (A16). This completes the proof. 



(A21) 



t( 4 )2_i_J 2 )5 
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We are now in a position to see why the algorithm given here works well. Lemma A.l 
and the setting-up (A5, A6), are put together to imply that we have sloved through 
Steps-1,11 the inverse problem 

H{4) ^ 4) ^) = Hi2) (-^^ (A22) 

for H( 2 \ where S £ image£>^(4) jr; (2) up to degree-4. This means that S(q,r]) coincides 
up to degree-4 with the generating function S(q, rj) for the inverse problem for G, so 
that the degree-4 inverse problem can be solved through Steps-1,11. 



